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Evolution of black-hole intermediate-mass X-ray binaries: 
the influence of a circumbinary disc 
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ABSTRACT 

[justham, Ra ppaport k, Podsiadlowsk] l)2006D recently suggested that black-hole low- 
mass X-ray binaries (BHLMXBs) with short orbital periods may have evolved from 
black- hole intermediate-mass X-ray binaries (BHIMXBs). In their model the secon- 
daries in BHIMXBs are assumed to possess anomalously high magnetic fields, so that 
magnetic braking can lead to substantial loss of angular momentum. In this paper we 
propose an alternative mechanism for orbital angular momentum loss in BHIMXBs. 
We assume that a small fraction 5 of the transferred mass from the donor star form 
a circumbinary disc surrounding the binary system. The tidal torques exerted by 
the disc can effectively drain orbital angular momentum from the binary. We have 
numerically calculated the evolutionary sequences of BHIMXBs, to examine the in- 
fluence of the circumbinary disc on the binary evolution. Our results indicate when 
6 < 0.01 — 0.1 (depending on the initial orbital periods), the circumbinary disc can 
cause secular orbital shrinking, leading to the formation of compact BHLMXBs, oth- 
erwise the orbits always expand during the evolution. This scenario also suggests the 
possible existence of luminous, persistent BHLMXBs, but it suffers the same problem 
as in Justham. Rappaport & Podsiadlowski (2006) that, the predicted effective tem- 
peratures of the donor stars are significantly higher than those of the observed donor 
stars in BHLMXBs. 
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1 INTRODUCTION 

If the gravitational mass of a compact star exceeds 
the maximum value ~ 2 — SM© of a neutron star 
(e.g. iRhoades fc Ruffinil Il974l) . this object should be 
taken as a black hole (BH) candidate. There ex- 
ist currently around twenty s t ellar-mass BH candidates 
(see iLee. Brown fc Wiier3l2002l: lOrosz et alJl2002l : ICasarej 
1200(11 . all of them are located in binary systems where their 
dynamical masses can be availably estimated. Nine of these 
systems are defined as compact BH X-ray binaries (BHXBs) 
with short orbital periods (< . 5 d) and low-mass donors 
(< IMrT.I jLee. Brown fc Wiiersll2002l: iRitter fc Kolbll2003t 
|Pods iadiaws ki^R^p^or^^laDn200^ It was estimated 
that there ma y exist > 1000 compact BHXBs in the Galaxy 
jWiicrs 1996: IRomaniill998ri . 

The short orbital periods of BH low-mass X-ray bina- 
ries (BHLMXBs) imply that they must have undergone sec- 
ular orbital angular momentum loss. If their progenitor sys- 
tems contains a low-mass secondary initially, it is not clear 
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whether the secondary star has enough energy to eject the 
envelope of the black hole progenitor during the common 
envelope evolution phase (Po dsiadlowski. RaDoaoort & Ha3 
l200,'j:l,Tustham. R.aPDaport fc Podsiadlowsk i 2006). This dif- 
ficulty can not be dealt with by assuming a intermediate- 
mass secondary, which have radiative envelope and are 
not expected to be subject to magnetic braking (|KawaleJ 
^88), since the binary orbits will be widen when mass 
is transferred from t he less mass iv e secondary to the 
more massive BH (see lHuandll96.'j^ . lEeeleton fc Verbuii^ 
(198(|) suggested that these systems may be evolved 
from triple sys tems via the merging of the two massive 
components. Pod siadlowski. Cannon fc R.eesI I0jj2,5) pro- 
posed that massive Thorne-Zytkow objects could lead 
to the birth of BHs with a low-mass donor star. 
I.Iustham. Rappaport fc Podsiadlowskil l)200fi|) summarized 
the previous proposals for the formation of BHLMXBs. They 
instead suggested that BHXBs containing intermediate- 
mass Ap and Bp donor stars, which possess strong magnetic 
fields, may be driven to compact BHLMXBs via magnetic 
braking with an irradiation-driven w ind from the donor star 
(see however. lYxmgelson et alJl200(J) . 

Are there additional angular momentum loss mecha- 
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nisms for BHXBs besides magnetic braking? Here we explore 
an alternative possibility. It is well known that some frac- 
tion of the transferred matter from the donor star may leave 
the s ystem in various ways during the mass exchange pro- 
cess ivan den Hewej^Q^. In partic ular, for cataclysmic 
variab les (CVs), ISpruit fc TaamI l^20nl^ and lXaam fc Spruiti 
i200ll) suggested that part of the outflow may be in the 
form of a circumbinary (CB) disc, and investigated the in- 
fluence of the CB disc on the evolution of CVs. The simi - 
lar idea was adopted bv lMever fc Mever-Hofmeisterl ll200lll 
for BHXBs to regulate the mass transfer rates, in which 
the CB disc could result from the remnants of the previous 
common-envelope phase. In our previous works we have in- 
vestigated the evolution of neutro n star (NS) LMXBs and 
Algol type binaries wit h a CB disc llChen. Li fc Wan j|2006t 
IChen. Li fc Qianll2006ll . These works indicate that the pres- 
ence of a CB disc can accelerate the evolution process of the 
binary systems, enhance the mass transfer rates, and lead 
to secular orbit shrinking under certain conditions. 

The aim of this work is to study the influence of the CB 
disc on the evolution of BH intermediate-mass X-ray bina- 
ries (BHIMXBs), and explore the possibility of BHIMXBs 
being the progenitor systems of short period BHLMXBs. 
The structure of this paper is as follows. In section 2 we 
describe the adopted orbital angular momentum loss mech- 
anisms in the evolution model of BHIMXBs. In section 3 we 
present the numerically calculated results for the evolution- 
ary sequences. We summarize and discuss the uncertainties 
in our model in section 4. 



2 DESCRIPTION OF THE MODEL 

Although s ome BHXBs may form dynamically in globu- 
lar clusters l|R.appaport et al]l200J^ . we only consider those 
evolved from the primordial massive binary systems. We fo- 
cus on the evolution of BHXBs consisting of an intermediate- 
mass donor star (of mass Aid ~ 3—5Mq) and a black hole (of 
mass Mbh). Neglecting the spin angular momenta of both 
components, we consider three types of angular momentum 
loss from the binary system (in all calculations, magnetic 
braking mechanism is not included), which are described as 
follows. 



2.1 Gravitational wave radiation 

For a compact BHXBs gravitational wave radiation is able 
to carry away orbital angular momentum effectively and lead 
to mass transfer. The angular m omentum loss rate due t o 
gravitational radiation is given bv llLandau fc LifshitJliol : 



32 G^/^ MjjiMjM^^^ 
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(1) 



2.2 Isotropic winds 

Mass and orbital angular momentum loss may occur during 
rapid mass transfer phase, since the mass accretion rate of 
a black hole is limited by Eddington mass-accretion rate^. 
For spherical accretion, this maximum accretion rate can 
be derived from the equation that gravity can balance the 
radiation pressure: 



MRdd = 



4t7GMbh 

KCTJE 



(2) 



where k — 0.2(1 + X)cm^g^^ is the electron scattering 
opacity for a composition with hydrogen mass fraction X 
(KioDcnhahn & Wciacrt 1990). The energy release efficiency 
77e of disc accretion onto black hole can be approximately 
written as: 
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as Mbh < VGMbh, where Mr^ is the initial m ass of the 
BH (see, e.g. lBardeenlll97ol: iKing fc Kolblll999^. Fr o m the 
above equations. IPodsiadlpwskL^ippaportfc HarJ (|200.'^) 

obtained the expression of the Eddington mass-accretion 
rate for a disc-fed BH: 
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(4) 

In our calculations we assume that the transferred mat- 
ter in excess of the Eddington accretion rate is ejected in the 
vicinity of the BH in the form of isotropic winds, carrying 
away the specific orbital angular momentum of the BH. We 
introduce the parameter / to describe the fraction of mass 
loss from the binary system, defined by the following rela- 
tions: 



and 



Mb 



M = Md/, 



(5) 



(6) 



where M is the total mass-loss rate from the system, Md 
is the mass transfer rate from the secondary, and Mbh is 
the mass accretion rate of the black hole, respectively. If 
the mass transfer rate \Md\ ^ -A^ecM, f ~ 0, else f — 1 + 
MEdd/Md. 

Based on the above equations and assumptions, the an- 
gular momentum loss rate via isotropic winds is 



where 



• ,MdMd - 
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J = a fiQ, 



(7) 



(8) 



is the total orbital angular momentum and = MdMsn/M 
is the reduced mass of the binary system, respectively. 



where G is the gravitational constant, c the speed of light, 
M = Mbh + Md the total mass of the binary, a the binary 
separation given by the Kepler's third law a = {GM /Q?)^^^ , 
where Q. is the orbital angular velocity of the binary system. 



^ We do not consider mass and angular momentum loss due to 
the wind mass loss from the donor star, which is not important 



for intermediate-mass stars. Irradiation effect on the stellar ^ 
is also not included. 
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2.3 CB disc 

We further assume that a smaU fraction 5 of the material 
overflowed from the donor star's Roche-lobe always feeds 
into the CB disc rather accretes onto the BH. Then Md in 
Eqs. (5)- (7) should be replaced by Md(l — 5), the BH mass 
accretion rate is Mbh = —(1 — ,f)(l — 5)Afd, and the total 
mass loss rate M = Md{f + S — fS). Tidal torques are then 
exerted on the CB disc extracti ng orbital angular mo men- 
tum from the binary system (sec Taam & Spruit 2001*). 

Similar as in JTaam & S pruit (2001), Spruit fc Taainl 
i200lfl . and lChen. Li fc QianI ll200d) . the viscous torque ex- 
erted at the inner edge ri of the CB disc^ on the BHXB can 
be shown to be : 

1/3 



Ti = JcB = jna^SAU 



(9) 



where 7^ = r^/a, t is the time since the onset of Roche lobe 
overflow (RLOF). In the s tandard a viscosity prescription 
ijShakura fc SunvaevlllQT^ . the viscous timescale t^i at the 
inner edge in the CB disc is given by: 

47^ 



where /3 = i/i/ri, a and Hi are the viscosity parameter and 
the scale height of disc, respectively. 

Combine Eqs. (9) and (10), the angular momentum loss 
rate via the CB disc can be written as: 

JcB = r?MdjcB, (11) 

where rj = S{t/tviY^^ , and job ~ iJ^/m is the specific or- 
bital angular momentum of the disc material at ri (e.g. 
fSob erman. Phinncy fc van den Hcuvcl 1997). 

The orbital evolution of BHXBs is governed by the 
change of the orbital angular momentum J of the system 
caused by the above mentioned three mechanisms. Differen- 
tiating Eq. (8) we get 
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where we have neglected the terms of fS. Neglect the angu- 
lar momentum loss due to gravitational wave radiation, the 
orbital evolution is governed by 



Jc: 



Md 
Md 
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(13) 



Insert Eq. (11) into Eq. (13), orbit shrinking, i.e. f2 > 0, will 
occur when 

Mbh 



V > 



Md _^2(/-K<S)Md 



3M 



In this paper we take ri/g 



(14) 
1.7 (see 



lArtvmowicz fc Lubowll994l : lMuno fc Mauerhanl200d) . then 
?7 > 0.77 when Md < Mbh, and /, 5 < 1 (see Fig. 1). 

In Fig. 1 we plot the expected relation between rj and 
Md when the orbital period is constant for BH masses 
A/bh ~ 7M0 , IOM0 , and 2OM0 , respectively. Generally a 
larger 77 is required for smaller Md and larger Mbh- This 
can be satisfied with an adequate mass input rate 5 and 
sufficiently long time t of RLOF. 

^ In this subsection we use the subscript i to denote quantities 
evaluated at the inner edge ri of the CB disc. 




[M J 

Figure 1. The critical rj = 5{t/t^{)^^^ as a function of the mass 
A/d of donor star for a constant binary orbit. The BH masses are 
taken to be 7 Mq (dashed curve), 10 Mq (solid curve), and 20 
A^Q (dotted curve), respectively. 



3 NUMERICAL RESULTS 

We have calculated the evolution of BHIMXBs adopt- 
ing an updated version o f the stellar evolution 
code developed bv [ Eg gletorJ lll97ll Fl OZ^ {see also 

3[ 



lHa,n. Podsia,dlow.ski fc EgdetoUll 9941: IPols etdiri99^. In 

the calculations we set initial solar chemical compositions 
{X = 0.7, Y = 0.28, and Z = 0.02) for the donor stars, and 
take the ratio of the mixing length to the pressure scale 
height to be 2.0. We include the afore-mentioned three 
types of orbital angular momentum loss mechanisms during 
the binary evolution. For the CB disc, we take a = 0.01 
and /3 = 0.03 (BcUc ct al. 2004). Setting the initial mass 
of the BH A/bh = WMq, we have performed the evolution 
calculations of BHXBs with an initial secondary of Md = 3, 
4, and 5 Mq. 

We first estimate the possible range of the fractional 
mass input rate into the CB disc, S, when the CB disc's irifiu- 
ence on the orbital evolution of BHIMXBs is effective. From 
Eq. (10) we find the viscous timescale t^i ~ 100 yr at the 
inner edge of the CB disc for a BH binary with an orbital pe- 
riod of ~ Id, and the duration of RLOF t ~ 10* yr, so we get 
S > Q.n{Ui/tY^^ ~ 0.008 from Eq. (14). For systems with 
an initial orbital period of ~ 3 d, t ~ 10^ yr, and 5 > 0.04. 
Obviously the longer t, the smaller 5 required. These values 
of 5 seem to be in agreement with |Taam fc Spruit (200^) 
who have shown that the CB disc is effective in draining 
orbital angular momentum from the system provided that 
5 exceeds about 0.01. In our calculations we find that CB 
discs with 5 >~ 0.01 — 0.1 (depending on the initial orbital 
periods) will cause runaway mass transfer with high mass 
transfer rates, while the evolutions with 5 significantly less 
than 0.01 is similar to t hose without CB disc s, as already 
noticed in CV evolution jTaam fc SpruidbOOlD . 

Figure 2 shows the examples of the evolutionary se- 
quence with an initial orbital period of 1 d. The solid and 
dotted lines in the figure correspond to two different values 
of (5 = (5.5 - 7) X 10"^ and 5 = (4 - 5) x 10"^ respectively. 
As can seen in the figure, mass transfer is initially driven by 
nuclear evolution of the secondary, leading to expansion of 
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Figure 2. Evolution of the donor masses Md (top), orbital periods Porb (middle), and mass transfer rates Md (bottom) for BHIMXBs 
with the initial orbital period Porb = 1 d. The left, middle, and right panels correspond to Md = SM© (solid curve 5 = 0.0055, dashed 
curve 5 = 0.004), AMq (solid curve <5 = 0.007, dashed curve 5 = 0.005), and 5M© (solid curve 5 = 0.007, dashed curve 5 = 0.005), 
respectively. Stable and unstable mass transfer processes are plotted with thick and thin curves, respectively. 




Figure 3. Evolution of the donor masses Md (left), orbital periods Porb (middle), and mass transfer rates Md (right) for BHIMXBs 
with an initial orbital period Porb = 3 d. The solid, dashed, and dotted curves correspond to the initial donor mass of SM© (<5 = 0.055), 
4Mq (5 = 0.07), and 5M0 (<5 = 0.09), respectively. 



the orbits. But this tendency is held up when the angular 
momentum loss via the CB disc becomes sufficiently strong. 
The mass transfer drops into a "plateau" phase at a rate 
~ lO~^M0yr~^ f or a few 10^ yr. These features are differ- 
ent from those in IPodsiadlowski. Rappaport fc Heinl (l2003l) 
for the standard evolution of BHIMXBs, in which the orbits 
always increase secularly. After that the mass transfer rates 
increase sharply as the secondary ascends the giant branch, 
but the final orbital evolution depends on the adopted value 
of 5. With the higher value of 5, a compact BHLMXB will 
be finally produced after a few 10* yr mass transfer. If mat- 



ter does not leave the CB disc, we can estimate from Fig. 2 
the mass of CB disc Mcb ~ 0.015, 0.028, O.OSSMq with the 
initial donor mass of 3, 4, 5Mq, respectively. 

The calculated results for slightly wider initial orbital 
period of 3 d are presented in Fig. 3. Generally the mass 
transfer rates are higher for higher donor masses and wider 
initial orbits - the peaks of the mass transfer rates can reach 
> lO~^M0yr~^. Within less than 1 Myr the donor masses 
decrease to be < IMq. Short orbital periods can also be 
attained, but generally a larger 5 (up to < 0.1) is required 
compared with those in Fig. 2. In this case the the CB disc 
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Figure 4. The evolution of tlie donor masses and orbital 
periods Porb for NSLMXBs with a initial donor mass of A/d = 
IMq, and an initial orbital period of Porb = 1 d. The solid and 
dashed curves correspond to 5 = 0.005 and 0.002, respectively. 
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to the thermal-viscous instability of accretion discs, when 
the mass transfer rates lie below a critical value Mcr so 
that the surface temperature at the outer edge of accre- 
tion disc is lower than the hydr ogen ionization tempera- 
ture jKing. Kolb fc Burderi[l993) . This critical mass trans- 
fer rate depends primarily on Mbh, Md, and Porb, and is 
given by (Ivan Paradiijll996l: lOubus et alj [l999): 

^""^■6^^° (lOM^j (iMej VTa) '"'^y' 

(15) 

In Fig. 2 the thick and thin curves correspond to stable 
and unstable mass transfer during the evolution of BHXBs. 
However, the mass transfer would be always stable for the 
wide orbit system with an initial orbital period of 3 d (see 
Fig. 3). It is interesting to see that the model BHLMXBs are 
likely to be persistent X-ray sources. This is not compati- 
ble with the observations of Galactic BHLMXBs, suggesting 
that S does not need to be constant but changes during the 
evolution (see discussion below) ^. 
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Figure 5. The evolutionary tracks of BHIMXBs in the effective 
temperature of the secondary vs. the orbital period plane when 
the initial orbital period Porb.i = 1 d. The solid, dashed, dotted, 
dash-dotted, and dash-dot-dotted curves correspond to the initial 
donor star masses of IMq (5 = 0.002), 1.5Mq (S = 0.0042), 
2Mq (S = 0.0048), 3Mq (S = 0.0055), and 4Mq (S = 0.007), 
respectively. The orbital periods of eight compact BHXBs and the 
earliest spectral types of the optical companions JCaaareEll200(j) 
are marked with crosses. 

is also more massive, Mcb ~ 0.15, 0.28, 0.45Mq with the 
initial donor mass of 3, 4, 5Mq, respectively. 

One may expect that the CB disc may play a similar 
role in the evolution of NSLMXBs. In Fig. 4, we plot the 
evolutionary sequence of NSLMXBs with a donor star of 
Aid = IMq and an initial orbital period of Porb = 1 d, 
to show the similar influence of changing S on the final or- 
bital evolution as in BH binaries. In our previous work on 
the formation of the bin ary radio pulsar PSR J1713-I-0747 
JChen. Li fc Wandl200d) we have found that for a given fi- 
nal orbital period, including the CB disc can allow a more 
evolved donor star at the onset of mass transfer, so that the 
X-ray hfetime can be considerably decreased. 

All of the BHLMXBs wi th an orbital period of < 0.5 
d should be transient sources JCasarejl2003l . It is currently 
believed that the transient behavior in X-ray binaries is due 



4 SUMMARY AND DISCUSSION 

In this paper we have examined the influence of a CB disc 
on the evolution of BHIMXBs. Assuming that a fraction S 
of the transferred material from the secondary forms a CB 
disc surrounding the binary system and extracts the orbital 
angular momentum from the binary, we have performed evo- 
lution calculations for BHXBs containing a donor star of 
mass in the range from 3 to 5 Mq with an initial orbital 
period of 1 d and 3 d. The calculations show that the evo- 
lution of BHIMXB is sensitive to 5. Generally a larger S 
is required for a wider initial orbit since the mass transfer 
itself during the evolution always causes orbital expansion. 
The results indicate that the orbits of BHXBs would show 
secular shrinkage when the values of 5 lie in the range of a 
few 10~^ — 10~^. Thus our CB disc scenario suggests a new 
evolutionary channel for the formation of BHLMXBs from 
BHIMXBs.' 

However, our results encounter the similar d ifficulty as 
in lJustham. Rappaport fc Podsiadlowskil i2006l) . that the 
calculated effective temperatures Tea are not consistent with 
those of the observed donor stars in BHLMXBs. We com- 
pare the calculated results of both BH-LMXBs and IMXBs 
with the observations in the Tcff - Porb diagram in Fig. 5. 
One can see that the observed results seem to be more con- 
sistent with the evolutionary tracks of origin al BHLMXBs. 
lJustham. Rappaport fc Podsiadlowskil ll200d) have already 
pointed out that this tempreature discrepancy seems to be 
a generic difficulty with any formation scenario that invokes 
primordially intermediate-mass donor stars. If this effective 
temperature problem can be solved, the CB disc mechanism 
may provide a plausible solution to the BHLMXB formation 



^ Chandra observations of the 12.6 hr ultraluminous X-ray source 
in the elliptical galaxy NGC3379 suggest th at the current on 
phase has lasted ~ 10 yr iFabbiano et al J2006l) . This source could 
be a soft X-ray transient with long duration of outbursts. How- 
ever, a persistent X-ray source as suggested by our calculations 
is an alternative possibility. 
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problem, without requiring anomalous magnetic fields in the 
donor stars. 

The mechanism feeding the CB disc is still unclear. It 
has been argued that during mass exchange in binary sys- 
tems, some of the lost matter which possesses high orbital 
angular momentum may form a disc surrou nding the binary 
syste m rather than leave the binary system ivan den Heuvej 
Il994) . This part of matter may come from the stellar wind 
from the donor star, wind and/or outflow from the accre- 
tion disc, or mass lost from the outer Lagrangian point. The 
values of 5 adopted here seem to rule out stellar wind as the 
origin of the CB disc, which requires a unreasonably high 
wind mass loss rate in intermediate- and low-mass stars. 
We speculate that the disc wind/outflow may play a more 
important role in feeding the CB disc, as the X-ray irradia- 
tion on the accretion disc in BHXBs may evaporate the disc 
much more efficiently than in CVs. This may explain why 
the values of 5 are about 2 — 3 orders of magnitu de larger 
tha n those for CV evolution l[S £ruit Taaml^2001^ . 

iDubus. Taam fc Sorxiitl (|2002h investigated the struc- 
ture and evolution of a geometrically thin CB disc to cal- 
culate its spectral energy distributions, and discussed the 
prospects for the detection of such d iscs in the i n frared 
and submillimeter wavelength regions. iDubus et alJ ll2004ll 
searched for excess mid-infrared emission due to CB disc ma- 
terial in CVs. But direct detection of the CB discs in CVs 
by infrared continuum studies has so far been elusive, partly 
because of the lack of accurate di sc atmosphere models. 
Recently, iMuno fc MaxierhanI l)200ft) studied the blackbody 
spectrum of BHLMXBs A 0620-00 and XTE J1118-f480, 
and found that the inferred excess mid-infrared emitting ar- 
eas are ~ 2 times larger than the binary orbital separations. 
Therefore, the detection of excess mid-infrared emission 
from these BHLMXBs provides evidence of the existence 
of CB disc around some BHLMXBs. These observations 
may set useful constraints on the evolution of the CB discs. 
The masses of the model CB disc are significantly larger 
than those (~ Mq) estimated bv iMuno fc MauerharJ 
i200 6f) . Combined with the non-detection of discs in CVs 
bylDubus ct al. (2004J, this fact suggests that the presence 
of the CB disc may not always accompany the RLOF pro- 
cesses. For example, (X-ray) nova bursts may destroy the 
CB disc before it has developed to be massive enough. Or 
perhaps 5 only needs to be high for a short time, while nor- 
mal magnetic braking takes over the CB disc during the 
majority of the evolution. The latter may partly account 
for the discrepancy between the lifetimes of BHLMXBs in 
Fig. 3 (~ 0.5 — 1.5 X 10^ yr) and those expected for observed 
BHLMXBs (~ 10® - 10® yr). Since the CB disc can promote 
mass transfer very efficiently, BHLMXBs are more likely to 
be observed when the value of 5 becomes extremely small. 
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